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CN ; Abstract 

A new fully kinetic system is proposed for modeling collisionless magnetic reconnection. The 
■ formulation relies on fundamental principles in Lagrangian dynamics, in which the electron inertia 

^< . is neglected in the expression of the Lagrangian, rather then enforcing a zero electron mass in the 

equations of motion. This is done upon splitting the electron velocity into its mean and fluctuating 
parts, so that the latter naturally produce the corresponding pressure tensor. If the electron heat 
" flux is neglected, then the strong electron magnetization limit yields a hybrid model for kinetic ions 

, and fluid electrons, in which the electron pressure tensor is frozen into the electron mean velocity. 

S: 

• Sh' Contents 

i-Sh . 1 Introduction 2 

i— i. 1.1 The kinetic model 2 

. 1.2 Electron pressure tensor dynamics 2 

^ ■ 1.3 Main content of the paper 3 

. 

^ " 2 Formulation of the model 3 

IT) " 

2.1 The action principle 4 

OO 1 

" 2.2 The variations 4 

<n : 

... 3 Electron pressure dynamics 5 

. ^ . 3.1 Hybrid model with frozen-in electron pressure 5 

^ . 3.2 Models with scalar electron pressure 7 

4 Stationary equations and constants of motion 7 

4.1 Stationary solutions and Harris' current sheets 7 

4.2 Total energy and constants of motion 8 

5 Conclusions 9 



1 



2 



1 Introduction 



A well known necessary condition for magnetic reconnection in plasmas is that the magnetic field 
lines are not frozen into the fluid flow. For example, while ideal MHD (with frozen-in magnetic field) 
fails to reproduce reconnection, a great advance is provided by resistive MHD, in which the frozen-in 
condition is broken by a finite resistivity. However, although reliable results are obtained by resistive 
MHD in several conditions, modeling collisionless reconnection requires extra effort, due to the special 
nature of the kinetic features underlying this phenomenon. Kinetic theory features first make their 
appearance in the dynamics of electrons, whose pressure tensor is often assumed to dominate over the 
inertia of their mean flow. In addition, the high energy levels of the ions lead to the necessity of a 
kinetic treatment also for these particles. 



1.1 The kinetic model 



In the light of the above arguments, fully kinetic simulations become necessary in modelling collision- 
less reconnection and they are based on a set of three equations, which may be conveniently written in 
terms of the ion distribution on phase space fi(x, v), the relative electron distribution / e (x, c) (where 
c = v — V e and V e is the electron mean velocity), and the magnetic induction field B. This paper 
proposes the following model for inertialess electrons: 



where 



dfe 

dt 
dfi 



+ (c + V e 



dfe 

d: 



+ 
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c x I — B - V x V ( 

m e . 



c-VV e + Vv 

m e n e 



- + v • — + — 

dt <9x rrii 



(V - V e ) X B + 



q e n e 



-V-P f 



dfi_ 
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Ve X B 
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fi d 3 v , q e n e V e = V x B - q { / v ft d 3 v 



dfe 

dc 







cc / d 3 c 



(1) 
(2) 
(3) 

(4) 



so that n e = J f e d 3 c is the electron density and P e is the electron pressure tensor. Here, the constants 
q e and q\ denote the electron and ion charge, respectively, and m e and rrii denote their corresponding 
masses. While the above dynamics of the magnetic field (3) and the ion probability density (2) is 
commonly found in the literature and it is easily derived by neglecting inertial terms in the equation 
for the electron mean velocity V e , the kinetic equation (1) for the relative electron density exhibits 
the new Coriolis-type force c x u e , where u e = V x V e is the mean electron vorticity. This type of 
non-inertial forcing appears typically in moving frames (recall that f e is the phase-space density in the 
frame moving with V e ), as shown in [17] for the case of general magnetized plasmas. However, this 
term is lost in common models, precisely in the step when one neglects electron inertia, after splitting 
electron mean and fluctuation velocities. This means that the asymptotic process leading to Hall 
MHD does not retain important features, such as non-inertial forces, when one aims to account for 
the fluctuation kinetics. Retaining these important features requires a different method. In particular, 
equations (1)— (4) are derived in this paper by applying Lagrangian methods: the same methods that 
yielded MHD in [7] and to Hall MHD in [10]. 



1.2 Electron pressure tensor dynamics 

In a series of papers (see e.g. [19, 22, 13, 14, 21] and subsequent papers by the same authors), 
considerable effort has been dedicated to the question of whether a moment truncation is possible, 
that would give a satisfactory description of electron kinetics. Such a moment truncation must account 
for the pressure tensor dynamics, so that the simplest possible truncation would rely on the assumption 
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of a negligible electron heat flux. Under this hypothesis, equation (1) eventually yields 

^ + (V e ■ V)P e + (V • V e )F e + ^(BxP e -P e xB)+F e .Vy e + (P e • vv e ) T 
at rn e \ J ' 

- (u e x P e - P e x cj e ) = . (5) 

V * ' 

Coriolis force terms 

Here, the new Coriolis term persists as it gives non-zero contribution in the pressure tensor dynamics. 
More importantly, one observes that in the strong electron magnetization limit BxP-PxB~0, 
the electron pressure tensor density is frozen into the electron mean velocity, that is 

j t (P ifc (qt,t) d</£ dq^ d 3 q t ) = along ^ = V e (q t , t) . 

This means that the electron velocity correlations are constant along the electron mean flow. Although, 
this is a nice picture, we should not forget that it assumes a negligible electron heat flux. The question 
whether this assumption is justifiable is still open, with results that basically depend on the particular 
situation that is being considered each time [14]. 

1.3 Main content of the paper 

1. Section 2 formulates the model by applying the hypothesis of negligible inertia for the electron 
mean flow. Instead of inserting this hypothesis in the equations of motion, this is done in the 
variational framework by using standard techniques in the theory of continuum systems with 
Lagrangian labels [7]. 

2. Section 3 analyzes the consequences of the model formulated in Section 2 in terms of the electron 
pressure tensor. It is shown how the strong electron magnetization limit yields the frozen- 
in condition for the electron pressure tensor, as long as the heat flux contribution may be 
neglected. Also, it is emphasized how this new frozen-in law for the electron pressure can be 
used to construct new hybrid models, which discard the information about higher order moments 
of the electron distribution. 

3. Section 4 presents the stationary equations in various cases and shows how the present model 
recovers the Harris current sheet solutions of the Maxwell- Vlasov system. In the same Section, 
a conserved total energy is presented explicitly, along with two families of constants of motion 
that are inherited from the Maxwell- Vlasov system. 

2 Formulation of the model 

The Lagrangian approach to continuum systems has a long standing tradition, whose most famous 
result is probably Arnold's formulation of Euler's fluid equation [1] in terms of geometry and symmetry. 
Later, this approach was extended to various compressible fluid models in different contexts [7]. In 
plasma kinetic theory, the Lagrangian approach appeared in the well known work by Low [12], which 
then was considered by Dewar in [4] . This approach achieved its most prominent result in Littlejohn's 
formulation of guiding center motion [11]. Later, this approach was pursued by many others [20, 2, 3]. 
When the variational approach arises from an action principle that is written in terms of purely 
Lagrangian labels, then this approach is often known as Euler-Poincare variational method [7]. For 
example, the variational formulations of ideal MHD [7] and Hall MHD [10] are exactly of this type. 
The present approach is mainly inspired by the results in [3, 10]. 
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2.1 The action principle 

In order to derive an appropriate action principle for modeling collisionless reconnection, let us write 
the Eulerian action functional as 



/ ^ ( y I 



+ ^jf e \u e + eV e \ 2 d 3 xd 3 c + q e j f e (u e + V e ) • Ad 3 cd 3 x - q e j v9/ e d 3 cd 3 x 
^j fi \ Ui -^ d 3 xd 3 v _H^j / e |5 e _ c | 2 d 3 xd 3 c-i f |Vx A| 2 d 3 x)dt = (6) 

where the notation is as in equations (l)-(4). In addition, (</?, A) denotes the electromagnetic poten- 
tials, e € [0,1] is a convenient parameter and the quantities 

Uj = Uj(x, v,t) , u e =u e (x,c,t) 

denote the position components of the Eulerian phase-space velocity for the two species, so that 
Xj = Uj(x, v, t) for the ions while x e = u e (x, c, t) + V e (x,t) for the electrons. In more generality, one 
defines the six-dimensional phase-space velocities Xj and X e so that, by a slight abuse of notation 
one may write 

(x, v)j = (uj,a,) =: Xj 

for the ions and 

(x,c) e = (u e + V e ,a e ) , with X e :=(u e ,a e ), 

for the electrons. Here, the accelerations aj and a e do not appear in the action principle (6), since the 
total energy cannot depend on the particle accelerations. Along the lines of [3], the terms in |uj — vj 2 
and |u e — c| 2 in (6) are used to constrain the velocities Uj and u e to the corresponding Eulerian 
coordinate so that 

u;(x,v,i)=v, u e (x, c,t) = c (7) 

Notice that expanding all terms containing Uj and u e in (6) yields Littlejohn's phase-space Lagrangian 
[ I ] for ions and electrons. This occurs because of the minus signs carried by the terms in — v| 2 
and |u e — c| 2 . These signs bring the present treatment very close to that in [3], although not identical. 

2.2 The variations 

At this point, one needs to take variations. Following [3] variations of the Lagrangian labels imply 
the following expressions for other variations 

5Xj = d t Y i + \Y i ,X i ], 5X e = d t Y e + [Y e , X e ] + [X Ve , Y e ] - [X w , X e ] (8) 
Sfi = -V • (fi Yi) , Sf e = -V • (f e Y e )-V- (f e X w ) (9) 

where [P, R] = (P • V)R — (R • V)P is the vector field commutator, while Y"i(x, v, t), Y e (x, c, t) and 
W(x, t) are arbitrary vector fields vanishing at t\ and ti. Here the vector field X^ e is constructed 
from V e as Xy e = (V e , c ■ V x V e ) and analogously for Xv^. Notice that, upon introducing the 
absolute velocity v := c + V e , we may write by a slight abuse of ntation 

(X, v) e = (Ue, S e ) + (V e , C • VV e ) = X e + X Ve 

The physical interpretation for the emergence of these expressions goes back to the fact that X e is 
a phase-space vector field that is expressed in the frame moving with V e ; see also [9], where the 
same method is applied to hybrid MHD models. Consequently, once variations (8) are taken in (6), 
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expanding the terms in Ssn and 5a e yields (7), while expanding the terms in Sui and Su e yields the 
expressions of the accelerations a, and a e . Eventually, one finds 



Xj = v . 



Ill; 



(v - V e ) x B + — V • 

q e n e 



(10) 



X e = f c , c x f ^B - V x V e ) - 2c • W e + — V • P e ) (11) 

V \ m e J m e n e J 

In the next stage, following [10], arbitrary variations (Sip, 5 A.) yield the first two equations in (4), 
i.e. Ampere's law and charge neutrality 

V x B = n e V e + f yfi d 3 v , and n e = - f ft d 3 v (12) 

On the other hand, variations of V e are given by 

sv e = d t w + [w, v e ] 

At this point, one makes the crucial assumption of inertialess electrons in the mean velocity equation: 
this step amounts exactly to letting 

e -> 0, 

in the action functional (6), so that the electron mean flow contributions to the kinetic energy is 
neglected, while retaining its current contribution in Ampere's law (12). The limit e — > in the action 
(6) yields the following expression for the electric field E 

I '8A \ 1 ~ 

E = - — + = -V e x B + V • P e , (13) 

V ot ) q e n e 

Ordinarily, a choice of gauge would be necessary to specify the dynamics of A. For example, in [10] the 
hydrodynamic gauge <p + V e ■ A = is chosen. However, taking the curl of the above equation yields 
(3). Notice that setting e = 1 in (6) amounts to including the electron inertial terms dtV e + (V e - V)V e 
in the right hand side of the second equality above. 

Finally, following [7, 3, 10], one recalls that (6) can be expressed in terms of Lagrangian labels and 
derives the kinetic equations 

Of- Of ~ ~ 

-jg + V • (fi Xi) = , -g- + V • (fe X e ) + V • (f e X V J = 

by taking the time derivative of the Lagrange-to-Euler map for the two probability densities. Then, 
using the expressions (10)-(11), one obtains the kinetic equations (1) and (2). 

3 Electron pressure dynamics 

3.1 Hybrid model with frozen-in electron pressure 

From the electron kinetics, one derives the evolution of the electron pressure tensor, which reads 

^ + (V e ■ V)P + (V • VeWe + ' W e + (P e " 

+ | — B - U e ) X P e - P e X { — H - U e ) = -V • Q e , 

where Q e = / ccc/ e d 3 c is the heat flux tensor. 

According to the above equation, the main problem concerning fully kinetic models for collision- 
less reconnection is that solving for two kinetic equations (ions and electrons) presents outstanding 
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computational difficulties [15] arising from the emergence of infinite moment hierarchies. This has led 
various authors to formulations of hybrid models, in which electron dynamics could be described by 
some kind of moment closure that retains the pressure tensor dynamics [5]. In particular, see also 
[19, 22, 13, 14, 21] and following papers by the same authors. In these works, the evolution of the 
electron pressure tensor dynamics is obtained by dropping the u e - and Q e -terms and upon replacing 
the B-terms as follows: 



+ (V e • V)P e + (V • V e 



+ 



vv e + 



•We) 



T 



IBI 



Pel 



where p e is the scalar electron pressure and I is the identity matrix, while r is a phenomenological 
time-scale introduced ad hoc to match the available data. Increasing values of r yield higher recon- 
nected flux, althuogh at some point singularities start to develop [22]. The above evolution equation 
yields reasonable results in many situations, although its underlying fundamental properties are not 
completely understood. Moreover, neglecting the heat flux tensor term may lead to physical inconsis- 
tencies [13], while a detailed study of its contribution is found in [14]. Still, most hybrid models for 
collisionless reconnection tend to disregard heat flux contributions. 

If we accept that the heat flux can be neglected V • Q e ~ 0, then, the model (l)-(3) (which retains 
the Coriolis force terms u e x P e — P e x cj e ) yields the electron pressure dynamics in section 1.2. 
Moreover, in the strong electron magnetization limit B x P e — P e x B ~ 0, the electron pressure tensor 
becomes frozen into the electron mean velocity, so that the electron velocity correlations are constant 
along the electron mean flow. When this happens, one can simply replace the kinetic equation (1) 
by (5). This yields a hybrid model in which ion kinetics (2) is coupled to an electron velocity V e 
transporting its own pressure tensor P e , along the lines of [19, 22, 13, 14, 21]. More explicitly, this 
procedure yields 



Of 



+ (V- V e )F e - [u e x 



dfi y dfj gj 
dt dx m,i 



(v-V e ) xB + — V 



X U e J + P e • We + (P e ■ = 

dfi 



dv 







<9B / ~ 
— = V X ( V e x B V • P £ 

dt \ q e n e 



(14) 
(15) 
(16) 



where 



q e n e = -q { / f { d 3 v , q e n e V e = VxB-g, v/, d 3 v , 



m e cc f d 3 c . 



(17) 



Remark 1 (Strong electron magnetization limit) The strong electron magnetization limit [21 J 

BxP e -P e xB~0 
is recovered in the present approach by neglecting the minimal coupling term 

q e J A • u e d 3 xd 3 c 

in the action (6), which then assumes that fluctuations are not directly coupled to the magnetic field. 

We point out that the Coriolis force terms cj e x Pe — Pe x u e are strictly necessary for a frozen-in 
electron pressure. Indeed, dropping these terms completely breaks this frozen-in condition. On the 
other hand, we should emphasize that this result is obtained by neglecting the heat flux contribution, 
which may not always be justifiable. The available data depend mainly on the particular situation 
that is being considered each time and no general statement is available in this regard. 



7 



3.2 Models with scalar electron pressure 

It is well known that a purely scalar electron pressure is not sufficient to generate reconnection. On 
the other hand, a scalar electron pressure is still used in many situations involving the generation of 
energetic particles. Examples are given by Field Reversed Configuration devices and certain space 
plasmas [16, 19]. In these contexts, hybrid models can be formulated, as described in [19]. 

When the kinetic features of electron dynamics are neglected and one is interested only in the 
quantities V e and n e , one only needs to replace the electron pressure tensor contribution V • P e in 
(2)-(3) by the pressure gradient Vp e , as it is given in terms of the electron density n e by an appropriate 
equation of state. The resulting model 



dfi y dfj qj 
dt 3x rrii 

dB 

~dt 



V e x B + 



1 



■v Pe 



<9v 







Vx(BxV f 



(18) 
(19) 



was also analyzed in [18] from the Hamiltonian viewpoint. Here, the scalar pressure arises from an 
internal energy function e(n e ) as p e := n^e'(n e ). 



4 Stationary equations and constants of motion 

In this section, we shall discuss a few properties of the stationary solutions of the model (l)-(3). While 
a nonlinear stability analysis requires knowledge of constants of motion that are not available a priori, 
some insight can be provided by the linearized equations. However, given the level of difficulty of 
the latter, this topic is left for further research in this direction. In this section, we shall make some 
comments on how the system (l)-(3) possesses Harris' current sheet solution and we shall show how 
the available constants of motion are insufficient for studying equilibria with non-zero current. 



4.1 Stationary solutions and Harris' current sheets 

Let us start our discussion by writing the stationary equations corresponding to (l)-(3). One has 



c + V e 



dfe 

<9x 



dfi . qi 

v ' + — 

ox rrii 

V x ( V P x B 



+ 



C X 



B — V X V e — c ■ VV. + 



-V-P, 



m P n P 



dfe 

dc 



(V - V e ) X B + 



1 



-V- 



q e n e 





d_h 
<9v 



with 



1 ' V' 



-qi Jfi d 3 v , 



q e n e V e = VxB-gJv/j d v . 



Now, following Harris' work [G], we look for stationary states such that E 
recalling (13), we have 



(20) 
(21) 
(22) 

(23) 

0. Then, upon 



V P . X B 



1 



q e n e 



-V-P, 



and we observe that the electron solution has the general form 



fe = fe[-\c + V, 



I 2 " VeP 

i 2 



f I I \ v \ 2 — - \v I 2 
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where we introduce v := c + V e . The above solution can be verified by replacing the above expression 
into equation (20), which may be written in terms of the canonical Poisson bracket {•, •} as 



1 



9fe 



/ e ,-|c|^+C- V e f + -^-( C + F e ) XB^ = 0, 



Notice that / e (x, v) and /j(x,v) satisfy 

9fe , 



v ■ — + I — v xB-v x V xV e + VV e ■ V e 

ox \m e 

dfi , qi <9/i n 

v • — H v x B • — = 0, 

ax m; av 



dfe 

dv 



(24) 
(25) 



where the terms in V e appearing in the electron equation are produced by the assumption of negligible 
inertia of the mean flow. 

Again, following [6], let us now restrict to the case when V e and Vi are spatially constant and 
V e = —Vi. Then, we obtain the usual Vlasov equations 



dfe . qe -r, dfe 

v ■ — 1 v x B • — — = 

ax m e Ov 

v — H v x B • — = 

ax 771; OV 



(26) 
(27) 



which, together with q e = e = —qi and 



n, 



/;d 3 v, 2n e V e = VxB, 



lead to the Harris current sheet solution in [6]. 



Remark 2 (Static equilibria) Notice that for static equilibria in which V e = f /;vd 3 v = 0, one 
has 



dfe 

ax 



Qe ^ 1 

— c x B + 

77i e q e n e 



dfi qi 
v ' IT H 

OX 777; 
V X ( — V - 5 



v x B + 



V 



-V • 







q e n e 



dfe 
dc 



av 







(28) 
(29) 
(30) 



and if we assume E = (q e n e ) 1 V • P e = 0, we have 

dfe , ? 



c — l^cxB.f = 
ox m e oc 

dfi , q% ^ „ dfi 

v ■ — H v x B • — = 

ax 777; av 



(31) 
(32) 



along with the zero-current relation V x B = 0. Here the equation (31) is different from (26), since 
the first involves only the fluctuation velocity. 



4.2 Total energy and constants of motion 

In the search for new physical models, it is essential that the energy is an exact constant of motion. For 
example, many models that are derived by making assumptions directly in the equations of motion do 
not possess this vital property (see [18] for a similar discussion on hybrid MHD models). In the case 
of system (l)-(3), exact conservation of energy is guaranteed by the underlying variational principle 
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(6). When this happens, the system is Hamiltonian, which means it conserves energy and it possesses 
a Poisson bracket structure. The latter property can be verified explicitly also in the present context, 
although this requires cumbersome calculations, which in turn do not add anything to the physical 
content of equations (l)-(3). Therefore the explicit identification of the Poisson bracket structure 
underlying (l)-(3) is left for future developments. Here, we shall focus only on the constants of 
motion. 

The following expression of the total energy remains constant under the dynamics (l)-(3): 

£( fl ,f e ,A) = ^ |/ J | v | 2 d 3 xd 3 v + ^ |/ e |c| 2 d 3 xd 3 c+i||VxA| 2 d 3 x, (33) 

as it can be easily checked by verifying that £ = explicitly. Notice that this is not the usual energy 
for the kinetic description of a hot plasma, mainly because only the fluctuation velocity c = v — V e 
enters in the expression of £. 

Other constants of motion are present in the dynamics (l)-(3). Actually, these are two separate 
families of constants that are defined as follows: 

a = /$i(/i)d 3 xd 3 v, C e = y$ e (/ e )d 3 xd 3 v 
where <l?j and $ e are arbitrary functions of their arguments. 

Remark 3 (Energy-Casimir method for nonlinear stability) Following [8], one can use the 

constants of motion to perform a nonlinear stability analysis for the system (l)-(3). In particular, 
this method identifies the equilibria by setting 

6(8 + d + C e ) = (34) 

and gives Lyapunov stability whenever 

5 2 {£ + d + C e ) > (35) 

However, the success of this method depends on the number of constants of motion that are available. It 
turns out that the three constants of motion found in this section restrict to consider only equilibria with 
zero current, thereby eliminating Harris ' solution from the treatment. Indeed, (34) yields V x B = 
at the equilibrium. At this point, one can try to find new constants of motion to enrich the stability 
analysis. On the other hand, this task can be difficult to reach and it requires finding the explicit 
Poisson bracket for the system (l)-(3). This will be subject of future work. 



5 Conclusions 

This paper has presented the new fully kinetic model (l)-(3) for collisionless reconnection, whose main 
novelty is the introduction of a Coriolis force term 

V7 T r df e 
- c x V x V e ■ — — 
oc 

in the electron kinetic equation. Unlike previous models, this term was obtained by neglecting the 
inertia of the electron mean flow in the variational principle underlying the Maxwell- Vlasov system. 
Perhaps non surprisingly, the Coriolis force arises from a change of reference, similarly to the results 
in [17]. 

The principal consequence of introducing the Coriolis force is that neglecting the heat flux contri- 
bution in the electron pressure tensor dynamics yields a frozen-in condition for the electron pressure 
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tensor itself, in presence of a strong electron magnetization. In turn, this frozen-in condition can be 
used to formulate new hybrid models, along the lines of [13] and subsequent papers (analogue hybrid 
models were also presented in [5]). The frozen-in condition in resulting hybrid models would confer the 
electron pressure an intrinsic Lagrangian meaning, which would open various possibilities concerning 
simulation techniques. 

The last part of the paper presented general considerations about the stationary equations for (1)- 
(3). Various specializations were considered and it was shown how the stationary equations comprise 
Harris' solution of the Maxwell- Vlasov system. Moreover, the explicit expression of the conserved 
total energy was provided, along with two more families of constants of motion. The question of how 
constants of motion can be used to perform a nonlinear stability analysis is left for future studies. This 
requires further constants that hopefully can be found from the explicit form of the Poisson bracket 
underlying (l)-(3). This Poisson bracket is the main object that underlies the Casimir method for 
Lyapunov stability and its identification can be be realized by cumbersome computations in the 
Hamiltonian framework. 
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